Thermal fluctuations of vortex clusters in quasi two-dimensional Bose-Einstein 

condensate 
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We study the thermal fluctuations of vortex positions in small vortex clusters in a harmonically 
trapped rotating Bose-Einstein condensate. It is shown that the order-disorder transition of two- 
shells clusters occurs via the decoupling of shells with respect to each other. The corresponding 
"melting" temperature depends stronly on the commensurability between numbers of vortices in 
shells. We show that "melting" can be achieved at experimentally attainable parameters and very 
low temperatures. Also studied is the effect of thermal fluctuations on vortices in an anisotropic 
trap with small quadrupole deformation. We show that thermal fluctuations lead to the decou- 
pling of a vortex cluster from the pinning potential produced by this deformation. The decoupling 
temperatures are estimated and strong commensurability effects are revealed. 

PACS numbers: 73.21. La, 73.30.+y, 79.60.Dp, 68.65.La, 36.40.Qv 



I. INTRODUCTION 

Properties of Bose-Einstein condensates (BEC) of al- 
kali atom gases attract a considerable current interest. 
Recent progress in this field allows for the creating of 
quasi-two-dimensional atomic gas either using ID opti- 
cal lattices or applying tight axial trapping 0, 0, S 0| • 
It is well-known that, according to the Mermin-Wagner- 
Hohcnberg theorem, Bose-Einstein condensation is im- 
possible in 2D homogeneous systems in the thermody- 
namic limit. However, a Bose-Einstein condensation at 
finite temperature becomes possible in a trapped gas. 

Recently, the Berezinskii-Kosterlitz-Thouless (BKT) 
transition associated with the creation of vortex- 
antivortexpairs was studied theoretically in 2D BEC 
clouds and it was shown that this transition can 

occur in the experimentally attainable range of parame- 
ters. For instance, according to Ref. [y], BKT transition 
can happen at T w 0.5T C for the number of particles 
N ~ 10 3 -j- 10 4 and realistic values of other parameters. 
These results demonstrate an importance of temperature 
effects in 2D BEC even at temperatures well below the 
critical one. At the same time, the effect of temperature 
on vortex lattices in BEC has not been studied yet, al- 
though fluctuation of positions of vortices should become 
considerable even at lower temperatures than those cor- 
responding to the BKT transition. Finally, experimental 
evidence for the BKT transition in trapped condensates 
was reported in Ref. Q- Recently, the effect of tem- 
perature on the vortex matter was analyzed in Ref. 
but in the strongly fluctuative regime at relatively high 
temperatures, when positions of vortices are random. 

It is well-known from the theory of superconductivity 
that thermal fluctuations can lead to the melting of flux 
line lattice. However, in real superconductors this usually 
happens only in the vicinity of the critical temperature. 
For the case of atomic BEC, critical temperature depends 
on the number of particles in the trap. Therefore, melting 
can occur at temperatures much lower than the critical 
one. In finite systems, fluctuations of vortex positions 



depend also on the number of vortices. In such systems, 
melting temperature is not a strictly defined quantity. 
In this characteristic temperature of the order- 

disorder transition ("melting") can be defined using the 
Lindemann criterion, see the discussion in |lCj . With in- 
creasing of the number of vortices, the fluctuations of the 
vortex positions are determined by elastic shear modulus 
of the system, i. e. by the Tkachenko modes studied 
in Ref. [TT| . However, when a vortex number is not 
large, quantization effects start to play a very impor- 
tant role and " melting" temperatures in this case can be 
much smaller than that for a larger system. Thermal 
fluctuations of the system of interacting point particles 
trapped by external potential were studied before in Rcfs. 
Oil H2, LL-E [HI mostly by using Monte-Carlo simula- 
tions. If there are not many particles (or vortices) in the 
system, in the ground state, they form a cluster consist- 
ing of shells. It was shown in Refs. [13, H3) LLil III! that 
with increasing the temperature, first, an order between 
different shells is destroyed and these shells become de- 
coupled with respect to each other. Only after this, with 
the sufficient increase of temperature, a radial disorder- 
ing of the cluster occurs. This leads to the hierarchy of 
"melting" temperatures, which depends dramatically on 
the symmetry of cluster and number of particles. 

In addition to thermal fluctuations, quantum ones can 
be significant in atomic condensate. In recent works 
[16l Il7j| , quantum and thermal fluctuations in finite vor- 
tex arrays in a one-dimensional optical lattice were con- 
sidered. See also Refs. 0, 0] for thermal fluctuations 
in spinor condensates. 

In the present paper, we study the intershell "melt- 
ing" of small vortex clusters in quasi-2Z? BEC at different 
number of vortices in the system. We consider the situ- 
ation, when a cluster consists of only two shells. First, 
we find the ground state configurations of vortices and 
then calculate deviations of vortices from their equilib- 
rium positions in a harmonic approximation. We show 
that, if the number of vortices in the inner and outer 
shell are not commensurate, deviations of shells with re- 
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spect to each other can be very significant even at low 
temperatures, T -C T c , and large number of particles in 
the system, and shells become decoupled with respect to 
each other thus leading to the disordering of the vortex 
cluster. Also studied is the role of thermal fluctuations 
on the small cluster, consisting on 2, 3, or 4 vortices, in 
a trap with a small quadrupole deformation, which acts 
as a source of orientational pinning for the cluster. 

The paper is organized as follows. In Section II we 
present our model, which allows one to find an energet- 
ically favorable vortex configuration in the 2D case and 
also to calculate semi-quantitatively deviations of vor- 
tex positions due to the thermal fluctuations. In Section 
III we study the intershell "melting" process in differ- 
ent two-shells vortex clusters and obtain a order-disorder 
transition temperature. In Section IV we analyze the ef- 
fect of thermal fluctuations on vortices in the trap with 
a small quadrupole deformation. We conclude in Sec. V. 



II. MODEL 



2D case: 



Hui± 
~~kT r 



02) 
N ' 



(3) 



where £(2) is a Riemann zeta function, yC(2) w 1.28. 
Eq. (3) remains accurate even for the case of interacting 
particles pcj . 



1. Ground state 

Now we present a method allowing one to find a ground 
state of the system, which corresponds to the certain vor- 
tex cluster, and deviations of vortices from their equilib- 
rium positions due to thermal fluctuations. 

In general case, ip can be represented as a Fourier 
expansion 



(4) 



Consider quasi-two-dimensional condensate with N 
particles confined by the radial harmonic trapping po- 
tential 



U(r) 



2 2 



(i) 



where uj± is a trapping frequency, m is the mass of the 
atom, and r is the radial coordinate. The system is ro- 
tated with the angular velocity to. In this paper, we 
restrict ourselves on the range of temperatures much 
smaller than T c . Therefore, we can neglect a noncon- 
densate contribution to the free energy of the system. 
Thus, the energy functional reads 

F = hu; ± N(T) J rdr J d^(||V^| 2 + y |</f + 2ng N \^\ 4 



■ i* 



(2) 



where the integration is performed over the area of the 
system, ip is the polar angle, N(T) is the number of con- 
densed atoms, u is the rotation frequency, qn = N \J^-^- 
is the interaction parameter, a and a z are the scattering 



length and oscillator length (a z — ■J ^-) in z-direction, 

which is kinematically frozen. Distances and rotation 
frequencies are measured in units of the radial oscilla- 
tor length and the trapping frequency, respectively. The 
normalization condition for the order parameter reads 
J rdr J dip \i[i\ 2 = 1. In this paper, we analyze the case 
of dilute BEC and take gN = 5, which corresponds to 
uj z /2it = 1.05 kHz at N = 1000 for 87 Rb (a ~ 5.3 
nm). Since we consider the range of low temperatures, 
T < 0.1T C , we can assume that N(T) ~ N. For the de- 
pendence of T c on N, we use the ideal gas result for the 



Let us denote a number of vortices in the system as u. If 
the superfluid phase in BEC has a g-fold symmetry, than 
only terms with Vs divisible by q survive in expansion 
(4). For instance, a vortex cluster consisting of a single 
ring of v vortices corresponds to the expansion (4) with 
I = 0, v, 2v, 3w, .... A two-shell cluster with v\ and vi 
vortices in the shells (wi+«2 = v), where V2 is divisible by 
v%, corresponds to the expansion (4) with I — 0, v\, 2v\, 
3i>i, ... If V2 is not divisible by V\, than, in general case, 
expansion (4) contains all harmonics. Typically, main 
contribution to the energy is given by just few harmonics 
and by taking into account approximately ten of them, 
one can find the energy of the system with very high ac- 
curacy provided that the number of vortices in the cloud 
is not too large, v < 10 -j- 20. 

In the limit of noninteracting gas (gN = 0), it follows 
from the Gross-Pitaevskii equation that each function fi 
coincides with the eigen function for the harmonic oscil- 
lator corresponding to the angular momentum I. These 

functions have the Gaussian profile ~ r l exp ^— . 

Therefore, one can assume that this Gaussian approxi- 
mation remains accurate in the case of weakly interacting 
dilute gas. The accuracy can be improved if we introduce 
a variational parameter Ri characterizing the spatial ex- 
tent of Finally, our ansatz for /; has a form: 



fl{r,Ci,R u 4> l )=C l 



exp 



2Rf 



(5) 



where C'i , Ri, and 4>i can be found from the condition of 
minimum of the energy (2), Ci is a real number. This ap- 
proach was used for the first time in Ref. |0J to evaluate 
energies and density plots of different vortex configura- 
tions. In Ref. 22], the simplified version of this method 
with fixed values of Ri = 1 was applied for the limit of 
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weakly interacting gas with taking into account up to 
9 terms in the expansion (4). In Ref. 0, the results 
for such approximate solutions to the Gross-Pitaevskii 
equation were compared with some known results of nu- 
merical solutions. A good accuracy of the ansatz was 
revealed. See also Ref. ,24| with the related approach. 
In Ref. [2{| , a version of this method was also used to cal- 
culate the energy of axially-symmetric vortex phases in 
spinor condensates with the comparison of the obtained 
results with numerical solutions, and a good agreement 
was found. Therefore, this method can be also applied 
for our problem and we expect that the results must be 
semi-quantitatively accurate and with the help of this 
model one can reveal the effect of symmetry of vortex 
cluster on melting temperatures and estimate values of 
those temperatures. 

Now we substitute Eqs. (4) and (5) to Eq. (2) and 
after the integration we obtain: 



I I l>k 

+4 ^ IlkkmClClC m Sl +mt 2k 
l>k>m 

x cos (cj>i + <j> m - 2<j>k) 

l>k>m>n 

X COS ((j>i +(f>k~ <Prn ~<t>n), (6) 



2. Thermal fluctuations: harmonic approximation 

After finding of ground state values of variational pa- 
rameters, one can calculate the equilibrium positions of 

vortices |t"o'\ </?o^}> 3 = 1) ■■• t S Dv numerical solution of 
equation 



(11) 



where we introduced a notation {p n } for the set of all 
variational parameters (Ri, Ci and (pi) and {fnj de- 
notes ground state values of these parameters. Fluctu- 
ations of p n , which are the degrees of freedom for the 
system in this model, lead to the fluctuations of vortex 
positions. We denote deviations of variational param- 
eters from their equilibrium values as 5p n and express 
deviations of vortices Sr^> and Sip^' through deviations 
of variational parameters in a linear approximation. Per- 
turbed positions of vortices are determined by the equa- 
tion 



5<pV\pV>+Sp n ) = 0. (12) 



Finally, deviation of the position of a given vortex from 
the equilibrium is 



A 



(J) 



where 



D(i) 



Sp r , 



(13) 



ai = -r(/ + 2) (1 + Rf) + nRfT(l + l)ul, 



*-lkmn 



2TT 2 g N T{ 



Rlkr. 



I + m + n + k 
2 

Rlkm 



Ri 
Rik in 



Rn 



R k 

Rlkran 



Rn 



DR. 



Ik it 



(7) 



(8) 



Sep 



(i) 



R 0) 



5Pn 



(14) 



Here and below repeated indices are summed; A 
Bn \ and are given by 



U) 

n j 



OPn 



Re 



dtp 

dPn 



Jin I 

or 



(15) 



Rlkr, 



V2{R~ 



Rl 



R 



-2\~2 



(9) 



T(l) is a gamma function. Normalization condition is 
now given by 



B<j)=Im(^-\Re (%£ 



dp r . 



dip 



Re 



dtp 

dp n 



(16) 



cfRtm +i) = i. 



(10) 



Values of parameters Ri, Ci and (pi can be found from 
the minimum of the energy (6) taking into account Eq. 
(10). For instance, for the axially-symmetric vortex-free 

state, C = 07^o> Rq = (1 + 2g) 1/4 , and d = at 
I 1. Note that the energy is proportional to frui±N at 
given values of w and gjy. 



D (j) = Im 



dip 
Or 



Re 



dip 
dp) 



Re 



dtp 
dr 



Im 



dtp 
dip 



(17) 

All the derivatives on p n , r, and (p in Eqs. (15)-(17) are 

taken at ground state values of parameters p n — pffi and 
space coordinates, corresponding to the equilibrium po- 



sition of a given vortex, r 



M 
'o ' 



The squared 



deviation of radial and polar coordinates of two vortices 
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labelled as j% and ji with respect to each other is given 
by 



G^SpmSpr, 



where 



A 



0* 1 ) 



A (J2 



DO 1 ) 




,4 



(i2) 



DO' 2 ) 



0i) 



DO' 1 ) DO' 2 




B 



02) 



DO 1 ) D0'2) 



(18) 



(19) 



(20) 



(21) 



In the same manner, we can express deviations of en- 
ergy from the ground state value as a quadratic function 
in terms of deviations of variational parameters: 



SF = E st Sp s Spt 



(22) 



where 



E, 



d 2 F 
dp s dpt ' 



(23) 



40) 



The derivatives here are also calculated at p n = pr, 

The averaged squared deviations of radial and polar 
coordinates of two vortices with respect to each other 
due to thermal fluctuations are given by 



6r (hh 



J d(5p)G < ^ l n 2) Sp m Spn exp(- yj:E st 5p s 5pt) 



J d(Sp) exp(~-^E st Sp s Sp. 



(24) 



__ J d(Sp) Jmn 2 $P m 5pn 6Xp(- -pp E st Sp s Sp t ) 

J d(Sp)exp(--p r E st Sp s 5p t ) 

(25) 

In general case, integrals in Eqs. (24) and (25) can not 
be calculated analytically, since matrix E st is not neces- 
sarily diagonal. Therefore, we have to switch to a new 
basis St = MSt, where M is a matrix, which diagonalizes 
quadratic form (22). Here and below we will use a matrix 
form for the equations. Quadratic forms (18), (19), (22) 
in the new basis can be written as 



5r thh 



(St) T pv^st, 



SF = {St) 1 QSt, 



(26) 



(27) 



(28) 



where Q = M T EM, = M T G^^M, and 

Rihh) = M T J^ lj ^M. Matrix Q must be diagonal and 
from this condition one can find matrix M numerically 
and then calculate pO'u'a) an d JjO'u'a), j n new basis, 
integrals in Eqs. (24) and (25) can be found analytically 
and finally we get: 



Sl Uth) 



= kT 



P, 



O1J2) 



Sif 



O1J2) 



kT- 



pOu'a) 



(29) 



(30) 



As usual in a harmonic approximation, average squares 
of deviations are proportional to the temperature. 

For the vortex cluster consisting of two shells we also 
introduce a quantity Atp, which has a sense of averaged 
displacement of vortex shells with respect to each other. 
It can be defined as a root of square displacement of pair 
of vortices from different shells averaged over all possible 
pairs of vortices: 



Aip - 



viv 2 



JX,J2 



1/2 



(31) 



Now, if we take into account Eq. (31) and the fact that 
energy in the ground state is proportional to Hu)±N, we 
obtain the following relation: 



t 1 ' 2 



d(g N ,u). 



(32) 



where t is a reduced temperature, t = T/T c ; function 
d(<7iv,u;) depends on interaction constant gjq and rota- 
tion speed u>. Of course, d(gpf,u>) is also very strongly 
dependent on the vortex cluster symmetry and in the 
next Section we will calculate it for some values of gjy 
and lo and vortex configurations. 

Note that the harmonic approximation remains accu- 
rate only if the deviations of the positions of vortices are 
much smaller than the characteristic distance between 
two neighboring vortices. The "melting" temperature 
can be defined through the Lindemann criterion. 



III. INTERSHELL MELTING OF VORTEX 
CLUSTERS 

If there are not many vortices in the system, they are 
situated in concentric shells. In the single vortex state, a 
vortex occupies the center of the cloud. If the number of 
vortices v is more than one, but less than six, vortices are 
arranged in one shell. With further increasing of number 
of vortices, one of the vortices jumps to the center of 
the cloud, whereas others are still situated in the single 
shell l2ol. Thermal fluctuations in these cases can lead 
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only to the radial displacements of vortex positions, since 
there is only one shell in the system. However, when the 
number of vortices is increased, they are arranged in two 
shells. For instance, it was shown in Ref. [21j, that in 
the phase with ten vortices, two of them are situated in 
the inner shell and eight are arranged in the outer shell. 
For phases with larger amount of vortices, their number 
in the inner shell can increase. 

Here, we consider the process of the intershell disor- 
dering in two-shells clusters containing 10, 11, 12, and 
13 vortices, respectively. It would be more convenient 
for the comparison to calculate "melting" temperatures 
for these configurations at the same value of rotation fre- 
quency. However, only one of these states can be a true 
ground state and if the system is not in a ground state, 
than sooner or later it will switch to the ground state 
due to thermal fluctuations. Therefore, we find "melt- 
ing" temperatures for different vortex configurations at 
different, but quite close to each other rotation speeds, 
which correspond to ground states of the given configu- 
ration. 

We choose the value of the gas parameter g^ = 5, as 
was explained in Section II, and find the ground states 
of the system. We have obtained that a two-shells vortex 
clusters consisting off = 10, 11, 12, and 13 vortices 
are energetically favorable in the vicinity of the point 
u) = 0.9. For instance, the ground state of the system 
is represented by phases with 10, 11, 12, and 13 vortices 
at u = 0.9, 0.91, 0.92, 0.94, respectively. In these cases, 
the inner shells contain V\ — 2, 3, 3, and 4 vortices, 
whereas the outer shells have v 2 = 8, 8, 9, and 9 vortices, 
respectively. Density plots for these vortex phases are 
shown in Fig. 1. Let's calculate deviations of vortex 
positions for these states. 

Using a technique presented in the previous Section, 
we found that if we increase the temperature from zero, 
at first, fluctuations of relative phases 4>i of different har- 
monics of the order parameter becomes important and 
deviations of the positions of vortices are almost entirely 
due to the fluctuations of <pi and not due to fluctuations 
of Ci and R\. This can be expected, since it is well- 
known that fluctuations of the phase of the order param- 
eter are more pronounced at relatively low temperatures 
and only at much higher temperatures an amplitude of 
the order parameter starts to fluctuate. Also, fluctu- 
ations of (f>i lead predominantly to azimuthal displace- 
ments of vortices, displacements in a radial direction are 
much smaller. This reflects the fact that the temperature 
of the intershell melting is much lower than that of the 
radial melting. We define an intershell melting tempera- 
ture t m eit of the cluster as a temperature, at which Aip is 
equal to 7^p, where 7 ~ 0.1 is a characteristic number 
from the Lindemann criterion: 

tmelt = ( ,f° , 7) VN- (33) 

\n 2 d{g N ,u) J 
Factor 360/«2 in Eq. (33) reflects the fact that the two- 



shell cluster is invariant under the rotation of shells with 
respect to each other on angle 360° /v 2 . We have cal- 
culated values of d for clusters with 10, 11, 12, and 13 
vortices. Our results are d{g N ,uj) w 608°, 3500°, 123°, 
and 810° for 10, 11, 12, and 13 vortex clusters, respec- 
tively. We can see that 10 and 12 vortex clusters are 
most stable among the analyzed configurations, and av- 
erage angle between the shells is less than in other cases. 
This is because the number of vortices in the outer shell 
v 2 is divisible by V\ . Intuitively, it is clear that the stabil- 
ity of a cluster with v 2 divisible by v\ depends also on the 
ratio v 2 /v\, since in the limit v 2 jv\ — > 1, each vortex in 
the inner shell corresponds to one vortex from the outer 
shell. Probably, this is a reason, why 12-vortex cluster is 
more stable than 10-vortex configuration. At the same 
time, 11 and 13- vortex clusters are the most unstable 
among considered here, since v 2 and v\ are incommensu- 
rate, and a deviation of shells with respect to each other 
is the largest. Note that with changing of oj with fixed 
qn, d(gN,uj) increases, in accordance with calculations 
[Tl| for Tkachenko modes. One can see from Eq. (34) 
and our estimates for d(gN,ui) that 12-vortex cluster is 
not melted and remains stable at N = 10 3 and t w 0.1, 
whereas in other cases a displacement angle between dif- 
ferent shells is comparable with the angle between the 
two neighboring vortices in the outer shell and therefore 
the shells are decoupled. The difference in melting tem- 
peratures for 12- and 11-vortex clusters is several orders 
of magnitude. 

Experimentally, melting of vortex clusters can be stud- 
ied by tuning of u> at fixed g^ and N. After obtain- 
ing a desirable vortex configuration, one can also tune 
T and reach a "melting" range of temperatures. Vortex 
positions can be found by the free expansion technique 
and after the repeating of this procedure one can deter- 
mine the average deviation of the vortex positions from 
the equilibrium. It is also possible to use a Bragg spec- 
troscopy for the systems containing much larger vortex 
arrays than those considered here. 



IV. VORTICES IN A TRAP WITH 
QUADRUPOLE DEFORMATION 

In this Section we consider the effect of thermal fluc- 
tuations on vortices in the trap with a quadrupole de- 
formation, which breaks a rotational symmetry. Such 
a deformation is often used in experiments to facilitate 
the creation of vortices. In fact, it introduces a prefer- 
able direction for the arrangements of vortices acting as a 
source of orientational pinning for a vortex cluster. The 
additional quadrupolar pinning potential is given by 

. . emeu 2 , r 2 cos2oj , , 

U quadr (r) = ^ £ (34) 

where e is a small coefficient, e<1. At zero tempera- 
ture quadrupole deformation of the trap potential fixes 
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azimuthal positions of vortices, whereas thermal fluctua- 
tions lead to displacements of vortices, which depend on 
e and T. Note that, in the case of a two-shell cluster at 
e = 0, considered in the previous Section, pinning cen- 
ters are created by each vortex shell for another shell, 
and the total rotational symmetry is preserved. Addi- 
tional potential (34) leads to the following contribution 
to the energy (6): 

P pm+4 Dm+6 

m V rn+2 ' JrL mJ 

xr(m + 3)cos(^ m+2 -<f> m ). (35) 

This term relates phases and amplitudes of different har- 
monics of the order parameter with angular momenta 
different by 2 from each other. Here we analyze the sit- 
uation, when there are v = 2, 3 or 4 vortices in the trap. 
Density plots for these vortex states are presented in Fig. 
2. The direction of minimum of trapping potential is ver- 
tical. We found that at small quadrupole deformation, 
e < 0.1, fluctuations of relative phases of the order pa- 
rameter Furrier harmonics are much stronger than fluc- 
tuations of their amplitudes, similarly to the case con- 
sidered in the previous Section. Therefore, we apply the 
same ideas to the present problem. Fluctuations of the 
relative phases lead mostly to the displacements of vor- 
tices in the azimuthal direction. These are so called scis- 
sors modes , w hich are responsible for such oscillations 

In the limit of small e, one can consider F quac i r as a per- 
turbation to the energy of the system with e = 0. If there 
are two vortices in the system, the nonperturbed order 
parameter contains all harmonics divisible by 2, and the 
main contribution to the energy is given by harmonics 
with I = and 2. Amplitudes of these harmonics are 
of e° order of magnitude, relative angles between them 
are fixed, and the energy is degenerate with respect to 
4>2, which reflects the fact that vortices can rotate freely 
(Goldstone mode). Quadrupole deformation connects 0o 
and <t>2, 4>2 and (f>4, 04 and <pe, etc.. As a result, F qua d r ~ 
e 1 and the angle of deviation of vortex cluster due to 
thermal fluctuations is given by 

t l/2 

A< ^ = £ r N l/A d ^N,^), (36) 

where d(gN,oj) is a function, independent on t; r = 0.5. 

If there are 4 vortices in the system, than at e = the 
order parameter contains all the harmonics divisible by 4. 
However, these harmonics cannot be related through the 
Eq. (35), since their angular momenta should differ by 2 
and not by 4. In this case, coupling of vortex cluster to 
the quadrupole deformations occurs in the next order of 
e. Namely, quadrupole deformation induces a harmonic 
with I — 2, whose amplitude is of the order of e 1 , and 
finally Fquadr ~ £ 2 i and we again arrive to Eq. (36) with 
t = 1. 



Now we consider the situation, when there are 3 vor- 
tices in the system. Again, a nonperturbed order pa- 
rameter consists of contributions with Vs divisible by 3, 
and their phases are not related by Eq. (35), as in the 
previous case with 4 vortices. In the next approximation 
with respect to e, quadrupole deformation induces other 
harmonics with all integer i's, and amplitudes of these 
harmonics are of the order of e . This is possible, be- 
cause any integer /, which is not divisible by 3, can be 
represented as 31 ± 2, and therefore it can be obtained 
by adding or substraction 2 from 31. However, it turns 
out that the energy even in the second order of e is again 
degenerate and this degeneracy is removed only in the 
next order of e. After all, F qua d r ~ £ 3 and Eqs. (36) is 
again valid but with r = 1.5. 

We see that the symmetry of the vortex configuration 
dictates the asymptotic behavior of the pinning energy 
and the average deviation of vortex cluster at e — > 0. 
The cluster with 2 vortices is most strongly pinned, 
whereas the cluster with 3 vortices is the most unsta- 
ble. This effect reflects the commcnsurability of the an- 
gular momenta of the quadrupole deformation (I = 2) 
and of the order parameter harmonics, responsible for 
vortices (I — 2, 3, 4). The strongest pinning is observed, 
when these momenta are equal to each other (two- vortex 
state), less stronger pinning, when these momenta are 
commensurate, but not equal (four- vortex phase), and 
the most weak pinning for incommensurate case (three- 
vortex state). Note that similar scaling relations for the 
frequences of scissors modes in 2 and 3 vortex states were 
obtained recently in Ref. [3(j. 

Next we calculate values of d(gN,^>) for 2, 3 and 4 
vortex cluster at gN = 5 and lu — 0.68, 0.75, and 0.78, 
respectively, where these vortex configurations are ener- 
getically favorable, according to our calculations. Our 
results are d{g N ,w) « 65°, 130°, and 11° for 2, 3, and 4 
vortex states, respectively. If the deviation angle Aip be- 
comes of the same order as the angle between two neigh- 
boring vortices in the cluster, 2tt/v, we will treat this 
vortex cluster as being depinned from the quadrupole 
deformation. By using Eq. (36) and this condition, one 
can easily obtain phase diagram in the (t, e) space. The 
example of the phase diagram is presented on Fig. 3 for 
= 5 and N — 1000. Each line determines the bound- 
ary between the pinned and unpinned vortex cluster for 
a given vortex configuration. Below these lines, vortex 
cluster is pinned and above it is unpinned. One can see 
that the region of stability of 2-vortex cluster is much 
broader than that for two other configurations. 



V. CONCLUSIONS 

In this paper, we studied the effect of thermal fluctu- 
ations on small vortex clusters in harmonically trapped 
rotating Bose-Einstein condensate at temperatures much 
lower than the critical temperature. First, we have con- 
sidered the clusters, consisting of two concentric shells 
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of vortices. These were 10, 11, 12, and 13 vortex struc- 
tures. We obtained that with increasing the temperature 
from zero, first, an order between the positions of vor- 
tices from different shells is destroyed, whereas the order 
within each shell is preserved. By using a Lindemann 
criterion, we defined the temperature, corresponding to 
the decoupling of two shells of vortices with respect to 
each other, which determines an order-disorder transi- 
tion. This "melting" temperature is strongly dependent 
on the commensurability of the number of vortices in 
shells; less commensurate clusters have lower "melting" 
temperature. For instance, the "melting" temperatures 
for the 11-vortex cluster consisting of two shells with 3 
and 8 vortices and for the commensurate 12-vortex clus- 
ter with 3 and 9 vortices in shells differ in several or- 
ders of magnitude. Intershell order-disorder transition 
can be observed at experimentally attainable range of 
parameters. We have shown that intershell melting in 
atomic condensates can occur at very low temperatures, 
especially for incommensurate clusters. 

Also studied are vortex clusters in the trap with small 
quadrupole deformation of the trapping potential, which 
acts as an orientational pinning center for vortices. We 



have analyzed the case of 2, 3 and 4 vortices in the sys- 
tem. We have demonstrate that the pinning energy de- 
pends very strongly on the number of vortices in the sys- 
tem. With tending the quadrupole deformation e to zero, 
the pinning energy becomes proportional to e v , where co- 
efficient v = 1, 3, and 2 for 2, 3, and 4 vortex configura- 
tions, respectively. This is due to the commensurability 
between the angular momenta / = 2 transferred to the 
system by the quadrupole deformation and I = 2, 3, and 
4, responsible for the creation of 2, 3, and 4 vortices, 
respectively. Average deviation angles between the vor- 
tex cluster and the trap anisotropy direction diverge in 
different power laws with tending e to zero for different 
vortex configurations, Aip ~ e T , with r = 0.5, 1.5, and 1 
for 2, 3, and 4 vortex clusters, respectively. 
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VI. FIGURE CAPTIONS 

Fig. 1 (Color online). Density plots for the states 
with 10, 11, 12, and 13 vortices. Dark spots correspond 
to vortices. 

Fig. 2 (Color online). Density plots for the states 



with 2,3, and 4 vortices in a trap with small quadrupolc 
deformation. Dark spots correspond to vortices. 

Fig. 3. The phase diagram of the 2, 3, and 4- vortex 
clusters in the trap with the quadrupole deformation at 
N = 1000. Above these lines the cluster is decoupled 
from the deformation and below it is coupled. 



